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ABSTRACT 

We present measurements of the specific star-formation rate (SSFR)-stellar mass relation 
for star-forming galaxies. Our deep spectroscopic samples are based on the Redshift One 
LDSS3 Emission line Survey, ROLES, and European Southern Observatory, ESO, public 
spectroscopy at z=l, and on the Sloan Digital Sky Survey (SDSS) at z—0.\. These datasets 
cover an equally deep mass range of 8.5<; log(Af,j/AfQ)<;ll at both epochs. We find that the 
SSFR-mass relation evolves in a way which is remarkably independent of stellar mass, as 
we previously found for the star-formation rate density (SFRD)-mass relation. However, we 
see a subtle upturn in SSFR-mass for the lowest mass galaxies (which may at least partly be 
driven by mass-incompleteness in the K-selected sample). This upturn is suggestive of greater 
evolution for lower mass galaxies, which may be explained by less massive galaxies forming 
their stars later and on longer timescales than higher mass galaxies, as implied by the 'cos- 
mic downsizing' scenario. Parameterising the e-folding timescale and formation redshift as 
simple functions of baryonic mass gives best fit parameterisations of t{AIij) cx A/^^^ "^ and 
1 + Zf{Mb) oc AI^^^^^. This subtle upturn is also seen in the SFR density (SFRD) as a function 
of stellar mass. At higher masses, such as those probed by previous surveys, the evolution in 
SSFR-mass is almost independent of stellar mass. At higher masses (log(Af*/Af0)>lO) the 
shapes of the cumulative cosmic SFRDs are very similar at both z=0.1 and z=1.0, both show- 
ing 70% of the total SFRD above a mass of log(Af,j/A/Q)>10. Mass functions are constructed 
for star-forming galaxies and found to evolve by only <35% between z=l and z=0.1 over 
the whole mass range. The evolution is such that the mass function decreases with increasing 
cosmic time, confirming that galaxies are leaving the star-forming sequence/blue cloud. The 
observational results are extended to z^2 by adding two recent Lyman break galaxy samples, 
and data at these three epochs (z^O.l, 1, 2) are compared with the GALFORM semi-analytic 
model of galaxy formation. GALFORM predicts an overall SFR density (SFRD) as a function 
of stellar mass in reasonable agreement with the observations. The star formation timescales 
inferred from 1/SSFR also give reasonable overall agreement, with the agreement becoming 
worse at the lowest and highest masses. The models do not reproduce the SSFR upturn seen in 
our data at low masses, where the effects of extinction and AGN feedback should be minimal 
and the comparison should be most robust. 
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1 INTRODUCTION 
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A useful metric for quantifying activity in a galaxy is the specific 
SFR (SSFR, SFR per unit stellar mass, SFR/M*). This represents 
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an efficiency of star-formation, since it measures the observed star- 
formation rate relative to tliat wliich it must have had in the past in 
order to build up the observed amount of stellar mass (e.g., Keimi- 
cutt et al. 1994; Brinchmann et al. 2004). In the local Universe, a 
correlation between the SSFR and stellar mass is observed (Brinch- 
mann et al. 2004), such that lower mass galaxies exhibit higher SS- 
FRs than their higher mass counterparts. A similar trend is also seen 
at higher redshifts (Brinchmann & EUis 2000; Feuhier et al. 2005; 
Elbaz et al. 2007; Noeske et al. 2007b; Zheng et al. 2007; Damen 
et al. 2009b), with the overall normalisation of this relation shifting 
to higher SSFRs at earlier cosmic times. One result of this evolv- 
ing relation is that the number of galaxies exceeding some SSFR 
threshold shifts from higher to lower mass galaxies with increasing 
cosmic time, generally referred to as 'cosmic downsizing' (Cowie 
et al. 1996), or 'downsizing in (S)SFR'\ 

The evolution of this relation between SSFR and mass rep- 
resents an important test of galaxy formation models, since both 
SFR and stellar mass depend critically on the prescriptions used 
for star-formation and mechanisms ("feedback") which act to sup- 
press it. Thus, observations which measure how SSFR and stellar 
mass are related and evolve with cosmic time are required to con- 
strain such models. Several observational limitations to building 
such a study exist. Possibly the most serious of these is the use of 
different SFR indicators which have different sensitivities to unob- 
scured (e.g., UV, [O II], Ha) and obscured star-formation (e.g., far 
infra-red emission, FIR); additional dependencies, such as metalUc- 
ity (e.g., [O II]); and contamination from non-starforming sources. 
The need to span a wide redshift range often means that different 
indicators must be used at different redshifts within a given survey, 
necessitating switching between different indicators, in which case 
systematic differences between the indicators might masquerade as 
evolution. 

Previous studies of SSFR-mass and its evolution have used ei- 
ther mixtures of different SFR indicators (e.g., Juneau et al. 2005; 
Noeske et al. 2007b), photometric redshifts which require SFR, 
redshift and dust extinction, etc. to be estimated, in an often de- 
generate way, from the same set of data (e.g., Feulner et al. 2005; 
Damen et al. 2009a; Zheng et al. 2007). Many studies also use 
stacking techniques, in which case contamination from a subset of 
objects (such as AGN) is often difficult to detect, and suffers from 
the limitation that only the average property of the stacked sample 
may be estimated (e.g., Zheng et al. 2007; Dunne et al. 2009; Pan- 
nella et al. 2009). Those which use spectroscopy to obtain redshifts 
(and/or spectroscopic SFRs) and consistent indicators at all red- 
shifts still only probe the most massive galaxies at higher redshifts 
(Elbaz et al. 2007; Cowie & Barger 2008; Maier et al. 2009) 

In Gilbank et al. (2010b, hereafter Paper 2) we presented a new 
survey (the Redshift One LDSS3 Emission lines Survey, ROLES) 
designed to target with mulit-object spectroscopy (MOS) galaxies 
an order-of-magnitude less massive than previously studied at z~l. 
We advocated the use of [O II] as a SFR indicator, empirically- 
corrected as a function of galaxy stellar mass in order to correct for 
extinction and other systematic effects with using [Oil]. This is 
useful as [O II] is more easily accessible to optical MOS instru- 
ments out to moderate redshifts z~ 1.2 and thus leads to rela- 



^ To distinguish it from other observations which point to this scenario 
from different avenues (such as the observation that locally-observed mas- 
sive galaxies exhibit older steUar populations than their less massive coun- 
terparts 'archaeological downsizing' (Kauffmami et al. 2003; Thomas et al. 
2005)). See Fontanot et al. (2009) for a comprehensive summary. 



tively efficient surveys. We showed by comparison with Balmer 
decrement-corrected Ha SFR that this mass-dependent correction 
is reliable locally (Gilbank et al. 2010a) and furthermore gives good 
agreement with the extinction estimated from IR measures. 

In this paper we construct a sample using a single SFR indi- 
cator from Q < z < 1^, spanning the widest possible mass range 
studied with spectroscopy (for redshifts and SFRs) over this full 
range. By combining our low mass spectroscopic survey at z=l 
(Paper 2) with higher mass data from the literature, and our lo- 
cal comparison sample built from SDSS (Gilbank et al. 2010a), 
we construct an unprecedented sample in this parameter space of 
mass-redshift-SI'R. In §2 we present the samples used; §3 presents 
our results for the SSFR-mass relation at z=0.1 and z=l, and for 
the mass functions of star-forming galaxies. §4 uses a toy model to 
describe the behaviour seen in the SSFR-mass relation, and com- 
pares several of our observational results with predictions from the 
GALFORM semi-analytic model of galaxy formation. A discussion 
of possible systematic errors is made in §5 and the conclusions are 
presented in §6. Throughout we assume a ACDM cosmology with 
0„ = 0.3, A = 0.7, and Ho = 70km s"^ Mpc"\ All magni- 
tudes use the AB system. All stellar masses and SFRs are calcu- 
lated assuming (or transforming to - see Appendix A) a Baldry & 
Glazebrook (2003, hereafter BG03) initial mass function (IMF). 



2 SAMPLE 
2.1 z=l data 

In order to minimise uncertainties caused by comparing different 
SFR indicators, a z~l sample is constructed from surveys using 
emission Une indicators for objects with spectroscopic redshifts. 
0.88 < z ^ 1.15 data are taken from the sample described in 
Paper 2. Of the two ROLES fields studied in paper 2 (the Faint 
Infra-Red Extragalactic Survey, FIRES, and the Chandra Deep 
Field South, CDFS) only the CDFS is considered here. This is 
due to the wealth of additional data available in the CDFS (in- 
cluding additional public spectroscopy covering higher mass galax- 
ies than probed by ROLES alone). In addition, the CDFS results 
dominate the ROLES statistics (199 low mass galaxies vs 86 for 
FIRES) due to the larger volume probed by the CDFS data. Briefly 
these data comprise multicolour photometry from FIREWORKS 
(Wuyts et al. 2008), and spectroscopic redshifts and [O II] SFRs 
from ROLES (Paper 2) for low mass (8.5< log(M,/M0)<9.5) 
galaxies. ROLES is a highly-complete, low-mass survey with a 
sampUng completeness >80% comprising 199 galaxies in the mass 
range 8.5< log(M,/Mo)<9.3 (Paper 2). The SFR is estimated 
from the Kennicutt (1998) relation (which assumes 1 mag of ex- 
tinction at Ha and an [O II]/Ha ratio of 0.5) converted to our BG03 
IMF. This is then corrected as a function of mass using the empir- 
ical mass-dependent correction of Gilbank et al. (2010a). Stellar 
masses are determined by SED-fitting of the photometry, at the 
spectroscopically-determined redshift, using a grid of PEGASE.2 
models and a BG03 IMF as described in Glazebrook et al. (2004). 

The high mass end in CDFS is supplemented with ESO pub- 
lic spectroscopy (Vanzella et al. 2008, hereafter referred to as the 
'F0RS2' data) for which we have measured [O II] fiuxes from their 

^ The results we present actually use Ha at z=0.1 to combat incomplete- 
ness for [O II] at high steUar masses in SDSS. For the low mass galaxies 
which are the primary aim of ROLES, our results are unchanged if we use 
[OH] directly instead of Ha at z=0.1 (c.f. Gilbank et al. 2010a). 
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Figure 1. (Vsos — *775) - *775 colour-magnitude diagram for galaxies 
with photometric redshifts at z~l (filled circles). Red open squares indicate 
galaxies with spectroscopic redshifts in the ROLES' redshift range from 
the FORS2 sample, and open blue diamonds indicate spectroscipically- 
confirmed z~ 1 galaxies from ROLES spectroscopy. The left panel shows a 
histogram of the colours, clearly indicating the bimodal nature of the galaxy 
colour distribution at z~l. The dashed horizontal line at Veoe — *775 = 1-3 
indicates the division between star-forming (blue cloud) and passive (red- 
sequence) galaxies. 



Figure 2. Mass vs observed A' -band magnitude for the z~l sample. Black 
points show all galaxies in the ROLES+FORS2 samples. In this work, 
only galaxies with masses down to log{Mt / Mq)=8.5 (vertical dashed red 
line) are considered and brighter than K=24.0 (horizontal dotted blue line). 
Galaxies with [O II]-SFR^;0.3Mq yr~^ are indicated by red diamonds. 
Galaxies with SED-fit SFRs^O.SMq yr~^ are indicated by blue squares. 
Solid black line shows running median of 15 galaxies in log{Mt / Mq)-K 
and dotted green line shows an extrapolation of this relation where the data 
are clearly complete. See text for discussion. 



ID spectra and computed masses using exactly the same FIRE- 
WORKS photometry as used for ROLES (Paper 2). The F0RS2 
sample is less complete than ROLES, comprising 73 galaxies in 
the mass range 9.7 < log(M*/M0)^11.2 and an average sam- 
pling completeness of ~30-40% (Paper 2) with a typical redshift 
success rate of 72% (Vanzella et al. 2008). 

We explore the effect of selecting star-forming galaxies by 
their blue colours. This is often done in spectroscopic samples to 
reject possible contamination from emission line objects where the 
emission is not from star-formation, such as LINERS (Yan et al. 
2006), and might be done in purely photometric surveys in or- 
der to reject non star-forming galaxies (see also §5.2). For these 
CDFS samples, the blue cloud is isolated by examining a colour 
bracketing the 4000A break for galaxies with photometric redshifts 
(photo-z's) from FIREWORKS within the ROLES redshift range 
of 0.88 < z ^ 1.15. Fig. 1 shows the (Veoe — *775) - 1775 
colour-magnitude diagram (CMD) for galaxies in the CDFS with 
mean photometric redshifts in the ROLES' redshift window (black 
points). The left panel shows a histogram of the (Veoe ^ 4775) 
colour which is clearly bimodal. The dashed horizontal line at 
(Vfjoe — 4775) ~ 1.3 indicates a dividing line between red and 
blue galaxies. The exact position of this line is not critical and 
any reasonable cut around this value leads to comparable results. 
Open red squares indicate galaxies with spectroscopic redshifts in 
our redshift window from the F0RS2 subsample. The majority of 
these objects inhabit the blue cloud, but there is also a significant 
subsample residing on the red-sequence. This latter category may 
comprise objects with very low levels of star-formation (since we 
have pre-selected objects with detectable [O II] emission), not de- 
tectable from broad-band colours, or possibly emission unrelated 
to star-formation such as LINERs (Yan et al. 2006). Rejecting red- 
sequence galaxies in this way reduces the sample size somewhat, 
to 60 rather than 73. For consistency, we recompute the complete- 
nesses in the same way as Paper 2, this time considering the blue 
sample separately and find that this makes negligible difference 



(except in the brightest K bin which was previously 70% com- 
plete), the overall completeness remaining at ~30-40% indepen- 
dent of K-band magnitude. 

Blue diamonds indicate z~l ROLES galaxies.'^ As can be 
seen, all but one object clearly belongs to the blue sequence (as 
might be expected given our selection criteria), so the ROLES sam- 
ple may effectively be considered as a blue cloud sample. Indeed, 
the lack of red objects in ROLES is strongly suggestive that these 
low mass galaxies possess little dust. 

2.1.1 Mass completeness ofz'^1 data 

Estimating the completeness in mass for this low mass galaxy sam- 
ple is difficult, since the spectroscopy is pushing the limit of the 
optical and NIR photometry. Our results probe down to masses of 
log(Af*/MQ)=8.5 at our magnitude limit of ^^=24.0. So, it is im- 
portant to check that K=24.0 is bright enough to sample the z~l 
star-forming galaxies of interest. In Paper 2, we used photomet- 
ric redshifts to show that galaxies potentially within our redshift 
range with SFRs below our [O II] selection threshold likely do not 
significantly contribute to the global SFRD at 1 . Pushing these 
photo-z's to fainter fluxes to attempt to locate galaxies below our 
^^=24.0 limit with stellar masses potentially within our mass win- 
dow is complicated by the rapidly increasing photometric errors 
below this limit. However, instead of resorting to the photo-z's we 
can use the spectroscopic results to search for potential selection 

It should be noted that some of the spectroscopic redshift slice-members 
(open symbols) do not correspond to galaxies with a best-fit photometric 
redshift in this slice (filled circles). The photo-z selected galaxies here only 
consider the value of the most-likely redshift, as we just require a repre- 
sentative z~l CMD to select red vs blue galaxies. In the analysis to assign 
spectroscopic redshifts, we consider the full photo-z PDFs (paper 2) and so 
galaxies with peak photo-z's outside this window may still have a signifi- 
cant probability of belonging to this slice. 
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biases within our sample. Fig. 2 sliows the observed /f-band mag- 
nitude as a function of fitted stellar mass for all objects in the ini- 
tial ROLES+FORS2 sample (i.e., prior to applying magnitude and 
mass cuts). The general trend of K versus stellar mass (and its scat- 
ter) can be seen. At masses above log(Af,/Af0)~9, the distribu- 
tion is well-separated from the K=24 magnitude limit (horizontal 
dotted line). Thus, a linear extrapolation of the (median-smoothed) 
relation around this mass towards the log(M,/MQ)=8.5 limit 
gives some indication of how the average relation should look in 
the absence of any possible selection bias. The average relation 
(measured by a running median of 15 objects) agrees reasonably 
well with the linear extrapolation from brighter magnitudes, indi- 
cating that the bulk of objects at the log(M,/M0)=8.5 mass limit 
are likely above the 7^=24 magnitude limit. Assuming the scatter 
in magnitude at a given mass remains constant, a visual estimate 
from this plot suggest that mass incompleteness at the magnitude 
Umit will be minimal. Attempting to quantify this, fading galax- 
ies at log(Af,/Af0)~9.5 to Iog(Af,/Af0)=8.5 assuming the same 
mass-to-light ratio (M/L) as the higher mass galaxies, a maximum 
of 25% of these objects would fall below the K-band magnitude 
cut-off. Thus we may estimate that we are at least 75% complete at 
the very lowest mass limit of ROLES. 

In order to check that our sample is not biased towards more 
actively star-forming galaxies near the magnitude limit of the sur- 
vey, the lowest [OII]-SFR galaxies (^O.SMq yr^^) are indicated 
in Fig. 2. It can be seen that these objects are distributed over a 
wide range in mass, and thus we do not see a 'pile up' of more- 
actively star-forming galaxies at the mass/magnitude limit of the 
survey. This tallies with the difference in observed i^-band mag- 
nitude at a given stellar mass from the extensive grid of PEGASE 
models used in the stellar mass fitting. At a given SFR for a model 
log(A//,/A^0)=8.5 galaxy (the mass limit of ROLES), for a wide 
range of metallicities, extinctions, and star-formation histories, the 
intrinsic scatter in if is « 0.5 mag or a factor of about 60% in 
M/L. The difference between a SFR of O.SMq yr"^ and O.IMq 
yr~^ (the former being the approximate ROLES limit and the latter 
thus being a factor of three below the nominal limit) corresponds 
to a systematic offset of 0.4 mag in observed K. Thus, while it is 
slightly easier to observe more actively star-forming galaxies in our 
sample, our empirical test shown in Fig. 2 suggests that we are not 
biased in this way. 

2.2 Local (z=0.1) data 

Local (z~0.1) data are taken from the SDSS Stripe 82 sample of 
Gilbank et al. (2010a). This sample was cut in redshift to 0.032 < 
z < 0.20 to ensure the inclusion of the [O II] line at the low red- 
shift end, and to avoid the effects of incompleteness/evolution at 
the high redshift end. For the present work we switch from using 
the [O II]-derived estimates of SFR to using the Balmer decrement- 
corrected Hq measurements directly. Gilbank et al. (2010a) showed 
that the [O II]-SFRs could be empirically corrected to agree statis- 
tically with the Ha-derived SFRs, but that the depth of the SDSS 
spectroscopy led to incompleteness for high mass galaxies when 
[O II] was used as the SFR indicator, due to the reduced sensitivity 
of [O II] to SFR with increasing M*. Since this work primarily re- 
quires high completeness, we choose to adopt Ha as our indicator 
at z~0, safe in the knowledge that it agrees, on average, with mass 
dependent empirically-corrected [O II]-SFR. To select galaxies be- 
longing to the blue cloud, we adopt the rest-frame {u — g) colour 
cut of Prescott et al. (2009) to separate red and blue galaxies, as 
used in Gilbank et al. (2010a). 



2.3 Specific star-formation rates 

As mentioned above, the SFR per unit stellar mass (specific SFR, 
SSFR) is a useful quantity for gauging the efficiency of star- 
formation. The inverse of SSFR defines a timescale for star- 
formation, i.e.. 



SSFR-^ octsFROC M^/Mi, 



(1) 



(where the SFR is the time derivative, A/, , of the stellar mass, Af« ), 
which is simply the time required for the galaxy to form its stel- 
lar mass, assuming its SFR remained constant. In the following, 
we will be considering the SSFR, and closely related quantities, of 
ensembles of galaxies. It is important to distinguish between the 
way the different quantities are calculated, depending on whether 
one is interested in the properties of the average galaxy, or the 
cosmic average of all galaxies. These distinctions are nicely de- 
scribed in Brinchmann et al. (2004) where they use r|^^ to denote 
the SSFR (rsFR) of the typical galaxy, and Tsfr to denote the 
volume-averaged equivalent. These can be written 

/SFR\ 



V 



-) = (- 



and 

V 
I'SFR 



PSFR 



(2) 



(3) 



Voliune-averaged quantities (Eqn. 3) are calculated following the 
method described in Paper 2. Briefly, this involves using the Vmax 

technique and weighting by sampling completeness estimated from 
detailed photometric redshift probability distribution functions. 



3 RESULTS 

3.1 The SSFR-mass relation 

SFRs and stellar masses from the z^O.l and z~l samples are com- 
bined to study the evolution of SSFR-mass in Fig. 3. The left 
panel shows this dataset without any colour pre-selection, whereas 
the right panel shows the same data after first removing the red- 
sequence galaxies, as described above. The SDSS (z~0.1) galaxies 
are shown as contours for clarity. Filled red circles show the mean 
in bins of stellar mass. The z^l data are shown as smaller filled 
circles, colour coded by dataset: black points are ROLES' galax- 
ies; smaller red filled circles show the data from the public FORS2 
spectroscopy in CDFS. Although individual data points are shown 
for each galaxy, the mean is calculated weighting each galaxy by 
a completeness and Vmax weight. Larger black circles show the 
median in bins of stellar mass for the z~l data. For the F0RS2 
data, the mean relation is not significantly lowered by including all 
galaxies with [O II] emission, rather than just blue cloud galaxies 
with [O II] emission. However, the colour cut makes a much larger 
difference to the mean SSFR measured in the SDSS. The effect of 
a colour cut to select star-forming galaxies will be discussed in §5. 
For now, we note that we prefer the sample removing red-sequence 
members (i.e., the right panel) and proceed to consider this. 

The dotted line shows the best-fit power law, log{SSFR) = 
— 0.42 log(Af*) — 5.49, which is consistent with all the points 
within their la errors, although the data are suggestive of a step- 
ping of the slope toward lower masses. The green lines indicate a 
toy 'tau' model which will be described in §4.1. 

In Fig. 3, we see that the log(SSFR)-log(mass) relation 
evolves to lower normalisation at lower redshifts (as found by many 
other works). At the high mass end, the relation evolves almost in 
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Figure 3. log(SSFR) versus log(stellar mass) for galaxies in our samples. The left panel shows data for galaxies of all colours (with significant emission 
lines) in the spectroscope samples; the right panel shows galaxies after application of the blue colour cut. Small, filled points show data at 0.88 < 2 ^ 1.15: 
black points are from ROLES, the higher mass data (smaller red filled points) are taken from ESO (FORS2) public spectroscopy. Contours denote SDSS 
data at z~0.1. Larger filled circles show mean SSFRs in bins of stellai' mass for the data in the two redshift ranges: blue symbols at z~l and red symbols 
at z~0.1. Dashed lines denote the three different SFRs, as annotated. In the right panel, the dotted black lines show best fit power laws to the z=l data, as 
described in the text. Green lines show best-fit tau models, discussed in §4.1, with parameters as indicated on plot. The green arrow indicates the possible 
effect of incompleteness, showing the systematic shift (head of arrow) caused by removing the 25% lowest SSFR points from the log(M»/Af0)~9.2 bin (tail 
of arrow). See text for details. 



parallel, as found by, e.g., Zheng et al. (2007). However, when the 
low mass ROLES data are included, an upturn in the relation at 
lower masses is seen. 

In order to verify that this upturn is significant, a boot- 
strap technique is used to perform 1000 realisations of the z~l 
data. The bootstrap error bars for each data point are compara- 
ble with, or smaller than, the \a sampling errors shown in the 
plot. Furthermore, if the approximate mid-points of the lower 
and higher mass samples, given by simply taking the centre-most 
bins of the ROLES and F0RS2 data (i.e., log(M*/M0)=9.O and 
log(M,/M0) = lO.6), then the lower mass sample is always sig- 
nificantly higher than than the higher mass bin in all but two of the 
1000 simulations. If a weighted mean of all three bins in the low 
and high mass samples is considered, then this would strengthen 
the significance further. The main uncertainty between the lower 
and higher mass data is likely to be systematic in nature (discussed 
in §5), and so we also consider the following bootstrap test. The 
null hypothesis that the data in the lowest mass bin is in fact the 
same as that of the highest mass ROLES bin, log(A/*/Af0)~9.2, 
(i.e., that the z=\ SSFR-mass relation is in fact flat and that 
the observed upturn is artificial) is adopted. Data points from the 
log(Af,/Af0)~9.2 bin are taken and bootstrap resampled, with 
the lowest 25% of these points removed each time (making the 
conservative assumption that the ~ 25% incompleteness, estimated 
in §2.1.1, systematically occurs for the lowest SSFR values). This 
generates a new mean SSFR value in each realisation which is ob- 
viously systematically higher than that without the lower points re- 
moved. In 1000 realisations, only seven are as high as the mean 
SSFR observed in the data for the lowest mass ROLES point. A 
limit to this bias may be estimated from the 95 percentile of the 
bootstrapped distribution. This is indicated by the green arrow in 
Fig. 3 and is approximately consistent with the lower la error 
bar of the lowest mass ROLES point. This shows that 5% of the 
time, the upturn seen in our data may be produced from a flat 
relation, with the lowest mass and lowest SFR galaxies dropping 
below the K-hanA limit of our sample. We emphasize again that 



this test is likely conservative, since it assumes that the distribu- 
tion of log(M*/Af0)~8.5 galaxies is exactly the same as those of 
log(M,/Af0)~9.2 galaxies, and that the 25% lowest SSFR galax- 
ies are removed in each bootstrap realisation. The fact that we mea- 
sure an upturn in this data is therefore intriguing but not cannot be 
completely dismissed as due to selection effects, and can only be 
addressed with even deeper K-band data in future mass-selected 
samples. 

The reason for this upturn will be discussed in §4.1. The low 
mass data should be the least uncertain in terms of total SFR esti- 
mated from [Oil], since these objects possess little extinction, as 
shown in §B2.2. Finally, it is worth noting that the combination 
of SDSS plus the ROLES+FORS2 data spans the largest range in 
mass and redshift probed by a uniformly-selected, spectroscopic 
sample. 



3.2 The mass function of star-forming galaxies 

The preceding has examined the relation between SSFR and stellar 
mass. It is instructive to look at how this mass evolves with time. 
Since ROLES is an emission line-selected survey, it is not possi- 
ble to measure the mass function for all galaxies at z= 1 from this 
sample. However, it is possible to measure the mass function of 
star-forming galaxies using the same Vmax method used to calcu- 
late the SFRD (Paper 2, eqn. 18), replacing the SFR of each galaxy 
with unity. The number density of blue, star-forming galaxies, (j), 
as a function of log(i\f./i\f0) is shown in Fig. 4. ROLES+FORS2 
CDFS data are shown as filled blue circles with error bars. Open red 
squares show the mass function of blue galaxies (0.75<z^l.OO) 
from Pozzetti et al. (2009) using zCOSMOS spectroscopic red- 
shifts. The left panel also shows the mass function for galaxies of 
all colours (0.75<z^l.00) from the zCOSMOS spectroscopic data 
of Pozzetti et al. (2009) (open squares). 

Pozzetti et al. (2009) classifiy star-forming galaxies by their 
photometric type from SED-fitting at their spectroscopic redshift 
(which is not quite the same as selecting star-forming galaxies by 
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□ □ all 0.75<z<1.00 (Pozzetti+ 10) 

a □ SFIng 0.75<z<1.00 (Pozzetti+ 10) 

• SFIng 0.88<z<1. 15 (R0LES+F0RS2) 
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Figure 4. Left panel: mass functions at z~l. Black diamonds and squares show data for galaxies of all types from Drory et al. (2005) and Pozzetti et al. 
(2009) respectively. Red squares denote (blue) star-forming galaxies from Pozzetti et al. (2009) and blue filled circles show (blue) star-forming galaxies from 
the CDFS from ROLES+FORS2 spectroscopy. Right panel: Comparison of the z~l mass functions with local values calculated from the Stripe 82 sample 
of Gilbank et al. (2010a). Black diamonds show the mass function for all galaxies, which may be compared with that computed from the NYU-VAGC in 
(Baldry et al. 2008). Blue and magenta lines indicate mass functions for star-forming galaxies (blue-selected) with SFRs>0.3 Mq yr^^ and >O.1M0 yr^^, 
respectively. Note that the depth of the ROLES data at z= 1 reaches the same mass limit as in the local data. 



their [O II] emission and blue colours, as we have done). The re- 
sults shown in Fig. 4 uses their SED-fit late type galaxies (SED- 
LTG) class (from their fig. 11), which should most closely resem- 
ble our blue star-forming selection. Our results are consistent with 
theirs, within the errors. Our points are systematically below theirs, 
but this small offset is consistent with cosmic variance expected 
from the CDFS volume (Paper 2), or this could represent differ- 
ences in the selection criteria. 

The right panel repeats a subsample of the z=l results and 
shows local comparison data. Open diamonds with error bars show 
the local galaxy mass function for galaxies of all types from the 
NYU-VAGC subsample of SDSS from Baldry et al. (2008). The 
black dashed histogram shows the equivalent mass function from 
the Stripe 82 sample (Gilbank et al. 2010a) computed for the 
present work. The good agreement between the two datasets taken 
from different subsamples of SDSS suggests that the Stripe 82 
dataset is complete in mass down to a limit of log(A'/«/MQ)~8.5 
for galaxies of all types (and likely higher for blue, star-forming 
galaxies which have lower M/L ratios than red sequence galaxies 
and thus are complete to lower masses in a given luminosity-limited 
sample). Hence the Stripe 82 data and ROLES may be fairly com- 
pared down to the ROLES mass limit of log(M,/M0)=8.5. The 
blue and magenta dashed histograms show the Stripe 82 data select- 
ing only blue cloud galaxies with SFRs>0. IMq yr~^ and O.3M0 
yr~^, respectively. 

The mass function of z~l star-forming galaxies from 
ROLES+FORS2 spectroscopy shows amazing agreement with the 
dashed blue histogram from SDSS local data. This histogram uses 
a SFR limit of 0. IM© yr~^ which is a factor of 3 lower than the 
O.3M0 yr~^ limit used at z~l. Arguably this is the fairest compar- 
ison since the average SFR density has decreased by a factor of ~3 
over this redshift range. This limit of 0. IMq yr~^ locally also cor- 
responds to the limited required for the SFRD to have converged 
(Gilbank et al. 2010a). In order to quantify the difference between 
the mass functions at the two epochs, we construct the cumulative 
sum of each, i.e., Sz = X/m-ii where Sz is the sum at redshift 
z of the individual galaxy's number densities, <j)m for galaxies more 
massive than Af, up to a maximum mass of log(Af,/AfQ) = l 1. 



This latter cut is placed so as to not be unduly affected by the 
small number of very massive galaxies in our sample. Examin- 
ing the ratio of these cumulative mass functions at the two epochs, 
S\.o/So.i shows that the z=1.0 star-forming mass function is ev- 
erywhere between ~6% and ~35% higher than the z=0. 1 star- 
forming mass function.'' Recently, Peng et al. (2010) presented an 
elegant, empirical picture of galaxy evolution, one of the key in- 
gredients of which is the apparent constancy of the mass function 
of star-forming galaxies over this redshift range and beyond. We 
have now measured this observationally over a much wider mass 
range than previously studied. The drop in the mass function of 
star-forming galaxies from z=1.0 to z=0. 1 confirms that galaxies 
must be leaving the star-forming sequence/blue cloud, perhaps at 
an even faster rate than that assumed by Peng et al. (2010). 



4 COMPARISON WITH MODELS 

In this section, the observed SSFR-mass relation is compared first 
with a simple toy model to parameterise galaxies' star formation 
histories as a function of their baryonic mass, and secondly with 
a recent state-of-the-art semi-analytical model of galaxy formation 
(GALFORM). The GALFORM predictions will also be confronted 
with additional observational results taken from Paper 2. 

4.1 Tau models 

One useful way to construct toy models for the star-formation his- 
tories of galaxies (e..g, Savaglio et al. 2005, Noeske et al. 2007a ) is 
to adopt a closed-box model and make the instantaneous recycling 
approximation. If we relate the SFR, ^t, to the mass of gas, Mg, via 
some star-formation efficiency factor, e, 

* = eMg, (4) 

* If instead of the O.IMq yr~^ SFR threshold (blue histogram) we adopt 
the 0.3Mq yr~^ SFR (magenta histogram), the difference grows steadily 
from Ri25% above log(Af,/A/Q)=9.5 to «65% at log{M, /Mq)=&.5. 
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and assume a fraction, R, of the stellar mass formed is instantly 
returned to the ISM, this leads to an exponentially declining SFR 
(e.g., Chariot & Bnizual 1991; Bruzual A. & Chariot 1993^) with 
an e-folding time of t where t = l/[e(l — 7?)], 



(5) 



if we write as a function of the baryonic mass, Mb, of the galaxy 
(initially all gaseous, i.e., Mh = Mg at « = Zf). is the formation 
redshift of the stars. 



t{z) - t{Zf ) 



(6) 



where t{z) is the cosmic time at which the galaxy is observed, and 
t{zf) is the cosmic time atZf. Then, 



'f{Mt,Zf) = eMt = 



Mt 



ril-R)- 
The above, together with the relation 

ft(z) 
Jt(zj) 

which gives 



M,{Mi,z) = 



R 



[1 - cxp(-r/r)], 



(7) 



(8) 



(9) 



specify M* and * in terms of the baryonic mass, e-folding time, 
and formation redshift of the galaxy's stellar population. R may be 
estimated from stellar population theory: R =0.56 for our BG03 
IMF (e.g., Hopkins & Beacom 2006). 

Noeske at al. (2007a) describe a model of 'staged galaxy for- 
mation' where they parameterise t and Zf d& power laws of the 
baryonic mass. 



(10) 



(a similar model was proposed by Savaglio et al. 2005 to explain 
the high redshift mass-metallicity relation), and 



l + Zf{Mb) = cpM^ 



(11) 



This methodology is applied to the data in Fig. 3 as a con- 
venient way to parameterise the behaviour of the star-formation 
properties of galaxies. Simultaneous fits to the mean values (large 
symbols) at the two epochs are performed using mpfit^ in IDL 
(Markwardt 2009), allowing the parameters a, Ca, P, and cp to 
float free. The best-fit parameters are indicated on Fig. 3, with the 
best fit curves overplotted as green solid lines for the z= 1 data and 
green dashed curves at z=0. 1 . The best fit exponents in the (pre- 
ferred) dataset from the right panel have values of a = —1.0 and 
P = 0.3, the same as those proposed by Noeske et al. (2007a), 
albeit with different normalisations in order to fit the lower over- 
all normalisation of the z~l data. Due to the high mass limit in 
the DEEP2/AEGIS data, Noeske et al. (2007a) could only mea- 
sure the average SSFR for the highest masses, (log (M, /M© ) ~ 1 1), 
where they were >95% complete. The low-mass end was an ap- 
proximate by-eye fit to the bulk of the data points. Interestingly, 
if we combine our highly complete (>80%) low-mass data with 
the DEEP2/ AEGIS data a very different dependence is found (See 
§B3). 



^ Such declining SFRs have also been historically referred to as "// mod- 
els", where /isFij = 1 — exp(— 1 Gyr/r) 



see http : //purl . com/net/mpf it 



With the values fitted in Fig. 3, the typical zj for a galaxy of 
log(M,/M0)=[9, 11] at z=l would be [1.7, 4.0] with r =[34, 4] 
Gyr; or for galaxies with these masses at z=0.1, the corresponding 
values would he Zf =[1.1, 3.7] and r =[64, 5] Gyr. 

So, the upturn seen in the Z'^l SSFR-mass relation towards 
lower masses can be explained by a non-zero /3, i.e., the forma- 
tion redshift was more recent for lower mass galaxies. Although a 
linear fit (P = 0) in log-log space is permitted (§3.1), the higher 
mass points all lie systematically above the linear fit (Fig. 3, right 
panel) and lower mass points systematically below, favouring the 
current parameterisation. With only the higher mass data used in 
other spectroscopic surveys, it is tmlikely that this upturn could be 
detected. 

Consider the implications for the mass function of star- 
forming galaxies (§3.2). In this toy model, individual star-forming 
galaxies are increasing in stellar mass: a typical log(M*/Mo)=9 
galaxy at z=l would increase its stellar mass by wO.5 dex by z=0.1 
and a Iog(A/,/A/,;;,) = l 1 galaxy at z=l by «0.1 dex. Some fraction 
of galaxies must also have their star-formation terminated and leave 
the blue cloud (Arnouts et al. 2007; Faber et al. 2007). These effects 
would act so as to shift the exponential cut-off of the mass function 
towards higher masses (and to shift the low mass tail toward higher 
masses by a larger amount), and to decrease the overall normalisa- 
tion of the mass function, respectively. Therefore processes must be 
occurring in order to balance these effects, preserving the apparent 
constancy of the mass function (Peng et al. 2010). 

4.2 The GALFORM model 

In order to gain some insight into the processes which might be re- 
sponsible for the above results, we compare the observations with 
a recent semi-analytic model of galaxy formation, the GALFORM 
model of Bower et al. (2006), implemented within the Millen- 
nium A'^-body simulation. This is a version of the Durham semi- 
analytic model which includes feedback from AGN to quench cool- 
ing within massive haloes. 

Firstly we add additional data in order to test the model over a 
wider redshift baseline (§4.2.1), and recast the SSFR-mass results 
in a maimer more directly related to the predictions made by the 
model (§4.2.2). 



4.2.1 Additional data, z=2 

To the z = 1 and z = 0.1 data shown in Fig. 5 (blue and black 
points, respectively), we add higher redshift data at z ~ 2 from 
two BX galaxy surveys (Sawicki & Thompson 2006; and Reddy 
& Steidel 2009, red squares and circles, respectively). These z ^ 2 
samples are used to extend the redshift baseline of our comparisons. 

It is important to note that the two stellar mass functions at 
z ~ 2 are derived in a different way from those at lower redshifts. 
The Sawicki (2010) MF is obtained from the BX galaxy LF using 
an empirical Mi7oo-stellar mass relation. This relation has been 
derived from a small but deep sample of « ~ 2 galaxies in the 
Hubble Deep Field for which stellar masses were estimated from 
multi-wavelength (B through H) SED fitting, and which takes the 
form M, = 0.68Mi7oo - 0.46 (Sawicki 2011). Sawicki (2011) 
use it to infer a stellar mass function from a revised version of the 
Sawicki & Thompson (2006) z ~ 2 UV LF that has been updated 
to take into account luminosity-dependent dust effects. We follow 
the same approach but apply the mass-luminosity relation to the 
original LF (Sawicki & Thompson 2006). This LF has a some- 
what shallower faint-end slope than the updated LF; together with 
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the Reddy & Steidel (2009) result (below), which comes from a 
steeply-rising z~2 LF, the two results span a conservatively wide 
plausible range of z~2 possibilities. At their faint/low-mass end, 
Reddy & Steidel (2009) follow a similar approach and an earlier 
version of the same HDF-based mass calibration. Specifically, at 
the low-luminosity end they correct their BX galaxy LF using their 
prescription for magnitude-dependent dust and then apply the Saw- 
icki et al. (2007; see also Sawicki 2011) SFR-stellar mass relation 
to arrive at the low-mass end of the MF; at the bright/massive end, 
they use stellar masses obtained from SED fitting for individual 
galaxies in their large sample of BX galaxies. 

The studies at z ~ 2 have several key differences from our 
study at lower redshift. The principal of these is that stellar masses 
are not derived for individual galaxies, as is the case at lower red- 
shifts, but rather stellar mass densities are calculated from UV LFs 
and empirical stellar mass - Mi 700 (or stellar-mass - SFR) relations 
— an approach that can propagate significant scatter (~ 0.3 dex) 
into the results. Additionally, these conversions are based on the 
assumption of a single, constant star formation history in the SED 
fitting. Finally, it should be kept in mind that both the z ^ 2 sam- 
ples rely on essentially the same stellar mass - A/1700 relation and 
the chief differences in the results stem from differences in their 
somewhat different UV luminosity functions. 

4.2.2 Timescale for growth 

As shown in Eqn. 1, the inverse of the SSFR defines a characteris- 
tic timescale. Similarly, the inverse of the volume-averaged SSFR 
(Eqn. 3) defines a timescale for the volume-averaged population 
(PMt / Psfr). Bower et al. (2006) studied the dimensionless ver- 
sion of this quantity which they denote as _R = pu^ / pSFRtniz), 
the ratio of the past to present SFR; where tniz) is the age of the 
Universe at redshift z. The GALFORM predictions from the Bower 
et al. (2006) model are plotted in Fig. 5. Our results are overplot- 
ted. These results are closely related to the inverse of SSFR plotted 
in Fig. 3. However there are some important differences. Note that 
since the Bower et al. (2006) prediction is for all galaxies rather 
than the star-forming subsample as we have defined it, our estimate 
of Pm^ / PSFR needs to include all galaxies. The change in psFR 
moving from a star-forming sample to galaxies of all types should 
be minimal. Locally, Salim et al. (2007) find that 95% of the SFRD 
is in blue galaxies (as we have defined our star-forming sequence), 
and we discussed in §3.1 the difference between using all z~l 
galaxies with non-zero SFRs versus those selected to be in the blue 
cloud. In order to recalculate pM^ for galaxies of all types, instead 
of using our measurement of the MF for only star-forming galaxies, 
we use the mass function of Drory et al. (2005). The difference be- 
tween this total mass function compared with the star-forming mass 
function is ~0.3 dex (although compatible within our broad uncer- 
tainties) at the high mass end, where we posses measurements of 
the SFR (log(Af*/MQ)<l 1), and smaller than this towards lower 
masses. The zCOSMOS-measured spectroscopic total mass func- 
tion (Pozzetti et al. 2009) agrees well with the Drory et al. (2005) 
measurement where they overlap at higher masses. Results at lower 
masses, where the photo-z estimate of the total mass function is not 
tested, should be regarded as the largest possible systematic error 
in the comparison with the models. 

Another, more minor, difference between Figs. 3 and 5 is the 
distinction between the quantities defined in Eqns. 2 and 3. For 
these data, this difference is smaller than the random errors in a 
given mass bin. Note that the values in Fig. 3 have been weighted 
by Vmax and the completeness of the surveys. 
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Figure 5. The timescale defined by Pm„ / PSFR as a function of stellar 
mass for all galaxies rather than just the star-forming sample. Red sym- 
bols show z~l data taken the SDSS Stripe 82 (Gilbank et al. 2010a); z~l 
data (blue symbols) are taken from ROLES (SFRD) with the stellar mass 
function from Droiy et al. (2005) for the same field; red symbols are z~2 
data taken from Sawicki & Thompson (2006) (open red circles, SFRD) and 
Reddy Sl Steidel (2009) (open red squares, SFRD) with the stellar mass 
functions taken from Droiy et al. (2005). Dotted horizontal lines show the 
Hubble time at the coiTesponding redshift, colour-coded as above. Filled 
lines are GALFORM model predictions. See text for details. 

This plot may be compared with the dimensionless ratio plot- 
ted in fig. 7 of Bower et al. (2006). The dimensionless quantities 
may simply be recovered by subtracting the log of the Hubble time, 
plotted as the dotted horizontal line, from the log of the timescale 
plotted as the curve. The main result of this comparison is the point 
at which galaxies of different masses make the transition from sig- 
nificant or increasing modes of star-formation to quiescent or de- 
clining SFRs. This is the mass at which the star-formation timescale 
at a given redshift equals the Hubble time at that redshift, i.e., the 
point at which the curve crosses the corresponding horizontal line. 
In the data, this transition mass occurs at log(Af*/AfQ)~l 1 at z=2 
but at log(A//./A//Q)«10.5 for z=1.0 and log(A/./AfQ)Ril0.2 at 
z=l. This is one manifestation of cosmic downsizing in the SSFR: 
more massive galaxies have transitioned earlier from active to qui- 
escent star-formation, relative to lower mass galaxies. This same 
trend is seen in the data. Indeed, the models and the data seem 
to agree best around the value of this transition mass. At z=0.1, 
at higher and lower masses, the model predicts somewhat higher 
values of the past-average to present SFR. We will return to this 
point below. Also, towards lower masses at z=1.0, the model is 
consistent with the individual measurement uncertainties'^ but lies 
systematically somewhat above the data. This means that the model 
does not reproduce the upturn at low masses in Fig. 3 which is driv- 
ing our fit of a Tau model with non-zero j3. At low masses, almost 
all galaxies should be star-forming at the level we can measure and 
so differences due to the ROLES data only selecting star-forming 
galaxies and GALFORM selecting all galaxies should be negligible. 
It is also worth noting that feedback from AGN is negligible in the 
models for low mass galaxies and so, whereas these effects could 
both play a part in the disagreement at the high mass end, the dis- 



EiTor bars are not plotted for clarity, but these may be approximately 
gauged from the mass function errors in Fig. 4 and the SFRD errors in 
Fig. 6. 
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agreement at the low mass end is pointing to inadequacies in the 
modelled physics for low mass galaxies. It is beyond the scope of 
this paper to examine in detail what processes might be responsi- 
ble, but we note that effects such as over-quenching for satellite 
galaxies (e.g., Gilbank & Balogh 2008) may be at least partly re- 
sponsible, and we point modellers to the current low mass data as 
an important test of galaxy formation theories. 

4.2.3 The SFRD as a function of stellar mass 

Another test for the GALFORM model is the evolution of the SFRD 
as a function of stellar mass. These were calculated for z~0. 1 
in Gilbank et al. (2010a) and z~l in Paper 2 using the standard 
lA'^max technique as described in those papers. To these, the data 
from Sawicki & Thompson (2006) and Reddy & Steidel (2009) are 
again added, converting to SFRD-mass as described above. Fig. 6 
(left panel) shows the SFRD, psFR as a function of stellar mass 
in three different redshift bins (left panel). The different coloured 
points refer to data at the different redshifts (black: z-^O.l, blue: 
z~l, red: z~2). Again, the lowest redshift data (black histogram) 
are taken from Stripe 82 data, where psFR is calculated as de- 
scribed in Gilbank et al. (2010a). The z~l psFR are taken from 
the (empirically-corrected) [O II] SFRD in Paper 2. The z~2 psFR 
data are taken from the LBG samples of Sawicki & Thompson 
(2006) (filled red circles) and Reddy & Steidel (2009) (open red 
squares). The overplotted coloured curves show the prediction of 
the GALFORM model. Errors are Poisson errors from the number of 
objects and do not include uncertainties associated with assuming 
different star-formation histories, etc. (but they do encapsulate the 
measured scatter in the adopted SFR-mass relation). 

The right panel shows the same information, but now plotted 
as the ratio between the SFRD at the redshift of interest and that at 
z=0.1, psfr[z — 0.1). Solid curves show the model predictions 
relative to the z=0.1 model prediction, and the dashed curves show 
the predictions dividing by the z=0. 1 data. 

The overall normalisation of the GALFORM model SFRD at all 
three epochs is reasonably well-matched to the observations. The 
best agreement occurs for the lowest masses at z=0.l and z=1.0 
where the ROLES data are best measured. A larger disagreement 
occurs around log(Af,/A^0)« 10-10.5, where the model overpre- 
dicts the data at both epochs. Interestingly this is around the same 
place as the best agreement occurs in the timescale plot (Fig. 5). 
Since the timescale is simply the stellar mass function divided by 
the SFRD, this implies a corresponding systematic difference in the 
GALFORM mass function. Comparing the GALFORM mass func- 
tion with that of the Baldry et al. (2008) total stellar mass func- 
tion (plotted in Fig. 4), the former lies systematically above the 
latter at the highest and lowest masses log(Af,/Af0)<lO.5 and 
log(M»/MQ)>11.0, but within the la errors of each point. In- 
terestingly, both the data and model showed a downsizing trend 
in the star-formation timescale (as described above). However, we 
claimed in Paper 2 that no downsizing trend in SFRD was seen in 
our data, and that the SFRD at z=0. 1 could be fit by a renormalised 
version of the local SFRD. The GALFORM predictions are in good 
agreement with our data, where the effect of downsizing should be 
most obvious. Our highest mass datapoint is only marginally dis- 
crepant (<2a). Differences are most easily seen in the ratio plots 
in the right hand panel of Fig. 6. On this plot error bars have been 
omitted for clarity, but may be gauged from the left panel. The ratio 
of model SFRDs at z=1.0 to z=0.1 (solid blue line) is very nearly 
independent of mass, and it is only at log(M«/MQ)>l 1.2 where a 
significant difference is seen. For the z=2 data, both observational 
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Figure 7. Fractional cumulative SFRD as a function of stellar mass for the 
two spectroscopic surveys of ROLES+FORS2 (blue dashed line) and SDSS 
Stripe 82 (black solid line). At high masses, the two distributions are re- 
mai'kably similar. More than 70% of the SFRD occurs in galaxies more 
massive than log(Af« /M0)>1O at both epochs. See text for discussion. 



results indicate the peak of the SFRD lies towards lower masses 
than that suggested by the GALFORM model. Indeed, examining the 
right panel of Fig. 6, the shape of the SFRD-mass relation seems 
to have changed remarkably little since z~0. 1. 

In order to more closely look for differences between the 
SFRD-mass distributions at the different epochs. Fig. 7 shows the 
fractional cumulative SFRD in galaxies more massive than a given 
stellar mass. The distributions are only plotted for the spectroscopic 
(z=0.1 and z=1.0) data, where the measurements are most secure 
and uniform. At high masses (log(Af,/AfQ)>10), the two distri- 
butions are very similar. Indeed, almost 70% of the total SFRD at 
each epoch occurs above this mass limit. At log(Aft/Af0)~9.5, 
the z=l distribution drops below that of the local value. This may 
be due to inadequate sampling of galaxies in this mass range, as 
this is the point around which the F0RS2 data end and the ROLES 
data begin. At the low mass limits of the data, the ROLES distri- 
bution begins to exceed the local measurement, and the gradient 
in the former is much steeper. This reflects the upturns we see in 
Figs. 3 & 4. If the deficit around log(Af,/Af0)~9.5 is due to a 
paucity of these mass galaxies in our sample, then the excess mea- 
sured here is likely even higher than 5%. This shows that high mass 
galaxies (log(A/*/A/0)>IO) appear to have evolved equally as a 
function of stellar mass between z=0. 1 and z=I.O, in SSFR (and 
its inverse, timescale) and SFRD, whereas low mass galaxies have 
evolved more (as shown by Figs. 3 and 7). The strength of any 
statement on this differential evolution is limited by the size of the 
current high mass sample. 



4.3 Other semi-analytic model comparisons in the literature 

Other works have compared their observations with different semi- 
analytic models applied to the Millennium simulation. Elbaz et al. 
(2007) used UV+24^m SFRs and spectroscopic redshifts in the 
GOODS-N and GOODS-S (CDFS) fields to test the models of Cro- 
ton et al. (2006). They found an increase in the normalisation of the 
SFR-M* relation between z~0 and z~ 1 of a factor of 6 in the data, 
and an increase of only a factor of 2 predicted by the models over 
the same range. Damen et al. (2009a) used UV+24/im data with 
FIREWORKS' photometric redshifts in the CDFS to compare with 
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Figure 6. Evolution of the SFRD as a function of stellar mass. Left panel shows observed SFRD at z~0. 1 (black histogram), z~l from [O II] (blue filled 
circles with error bars), and z~2 (red filled circles (Sawicki & Thompson 2006), open red squares (Reddy & Steidel 2009)). Solid curves are model predictions 
from galform where the redshifts are indicated by the same colours as the data. Right hand panel shows the same information, plotted as the ratio of the SFRD 
at the redshift of interest divided by the z=0. 1 SFRD. Solid curves show the model predictions relative to the z=0. 1 model prediction, and the dashed curves 
show the predictions dividing by the z=0. 1 data. Towards lower masses, the steeper gradient in the ROLES data reflects the low mass upturn seen in earlier 
plots. See text for details. 



the models of De Lucia & Blaizot (2007). They found that the log- 
arithmic increase of the SSFR with redshift is nearly independent 
of mass, for the higher mass bins which could be tested (down to a 
limit of log(A/*/Af0)~9.5 at z=0.9)*. Over the same mass range 
in the SDSS and FORS2 data (Fig. 3), our data appear consistent 
with the same result. They found that the same trend was seen in 
the De Lucia & Blaizot (2007) models, as it also appears to be in 
GALFORM from the parallel nature of the curves in Fig. 5. 

Damen et al. (2009a) also used the dimensionless growth rate 
{GRsf = (*/A/*)f//(2)), i.e., closely related to the inverse of the 
star formation timescale divided by the age of the Universe at the 
redshift of interest (and the inverse of Bower et al. 2006's R), as 
the observational quantity to compare with semi-analytic models 
of galaxy formation. The predictions were made by Guo & White 
(2008) which are based on the model of De Lucia & Blaizot (2007). 
As with the Bower et al. (2006) model, these semi-analytic pre- 
scriptions are applied to the Millennium simulation. Damen et al. 
(2009a) found that the growth rate was in good agreement with 
the models at z~0, but that the growth rate in the observations 
increased much faster than those of the models. The De Lucia & 
Blaizot (2007) model seems to predict a much flatter evolution in 
the growth rate/timescale than the Bower et al. (2006) model over 
the redshift range considered here (compare fig. 1 of Guo & White 
2008 with fig. 7 of Bower et al. 2006 model). Thus, although both 
semi-analytic models qualitatively reproduce the downsizing trend 
in cosmic star-formation, they both fail to reproduce important de- 
tails: the De Lucia & Blaizot (2007) predicts too slow a rate of 
the growth of stellar mass due to star-formation, although the val- 
ues for different galaxy masses at z=0 are quite accurately pre- 
dicted; whereas the (Bower et al. 2006) model predicts approxi- 
mately the correct growth rate for all but the most massive galaxies 
(log(M*/M0)>ll), underpredicting the local SFRD in high mass 

* Note that although we use the same FIREWORKS data, we do not re- 
quire such a high significance limit on the K-hand magnitude as Damen 
et al. (2009a), so we push the available photometry deeper by requiring an 
emission line to obtain a redshift, rather than fitting high-precision photom- 
etry to obtain a photometric redshift. 



galaxies. Conceivably these shortcomings could be due to simi- 
lar factors, however a detailed investigation into the semi-analytic 
models is beyond the scope of this paper. 

Damen et al. (2009a) and others (e.g., Zheng et al. (2007)) 
have claimed that the evolution in growth rate, or equivalently 
SSFR, is the same for galaxies of all masses. If our limiting mass 
was higher (more like log(M, /Mo)>9.5) we would have reached 
similar results. Our lower mass sample identifies an upturn in av- 
erage SSFR at lower masses which we interpret as evidence for a 
later formation epoch for lower mass galaxies in our toy-modelling. 



5 POSSIBLE SYSTEMATIC ERRORS 
5.1 SFR calibration 

Our measurements are most robust for low mass galaxies, for the 
reasons mentioned above. More uncertainty exists when results for 
high mass galaxies are included. The high mass data used comes 
from public spectroscopy of the CDFS. This is a relatively small 
field (which is especially problematic for the rarer, high mass 
galaxies) and these data have a sampling completeness of ~50% 
(and an overall completeness, including the redshift success rate, 
of «30%). However, the agreement between the CDFS data and 
the larger zCOSMOS survey converted to the same empirically- 
corrected [O II]-SFR tracer gives reassurance that cosmic vari- 
ance does not have a significant impact on these results (Appendix 
B3). Another concern is how well the empirically-corrected [O II]- 
luminosity traces the total (i.e. extinction-corrected) SFR. The cal- 
ibration is known to give good agreement between extinction- 
corrected Ha-SFRs at z~0.1 (Gilbank et al. 2010a), where the ex- 
tinction has been calibrated either by the Balmer decrement or IR 
luminosity. The disagreement between the DEEP2/ AEGIS SSFRs 
and the current dataset (see Appendix B3) at the high mass end 
might suggest that the empirically-corrected [O II]-SFRs are miss- 
ing extinguished star-formation detected in the 24/im observations 
(but note that the dicrepancy still persists for [O II] -only SFRs); or 
recent suggestions that the rest-frame \2^m. data overestimates the 
true SFR may be responsible (see Appendix B3). 
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The validity of a mass-dependent correction for [O II]-SFRs 
must be tested with an independent SFR indicator less-sensitive to 
extinction. It is an important open question (Kennicutt et al. 2009) 
whether extinction is more physically related to a galaxy's stellar 
mass or SFR (e.g., Hopkins et al. 2001; Brinchmann et al. 2004; 
StrazzuUo et al. 2010). Recently Garn & Best (2010) have pre- 
sented compelling evidence from a large SDSS sample that sug- 
gests stellar mass is the more fundamental, supporting our mass- 
dependent empirical correction for the effects of dust. In addition, 
[O II] luminosity is also sensitive to the effect of metallicity. As 
we discussed in paper 2, it is an open question how the effects of 
evolution in metallicty and dust may change our mass-dependent 
correction to the [O II] -SFR between that calibrated locally and 
that at higher redshift. It is possible that the evolutionary effects 
may in fact cancel out, since at z~l galactic metallicity is lower 
than locally (so [O II] luminosity increases for a given SFR), and 
dust extinction is likely higher (Tresse et al. 2007; Villar et al. 2008) 
(so [O II] luminosity decreases for a given SFR). Indeed, Mannucci 
et al. (2010) recently proposed a fundamental relation between stel- 
lar mass, metallicity and S¥K which is observed to remain con- 
stant, due to a conspiracy in the evolution of these parameters, up 
to z~2.5. We estimate the likely size of these effects below. 

For low mass galaxies, the situation is somewhat simplified as 

a) the dust extinction is extremely low, as we have shown above; 

b) theoretical models for the dependence of [O II] luminosity on 
SFR as a function of metallicity (Kewley & Dopita 2002) give 
good agreement with empirical calibrations (Gilbank et al. 2010a). 
Eqn. 15 from Kewley & Dopita (2002) gives the correction to SFR 
which must be applied to [O II] luminosity as a function of metal- 
licity. Taking the mass-metalUcity (M-Z) relations from Savaglio 
et al. (2005) for z~0.1 and z~0.7 (their eqns. 8 and 7 respec- 
tively), one can estimate the change in this theoretical correction 
between z~0.7 and z^O. 1 . For a log (M, / Mg ) ~9.0 galaxy (which 
is well sampled by the Savaglio et al. 2005 data at z~0.7), the 
evolution in this correction due to the evolution in metallicity is 
«5%. This is not large enough to explain the observed upturn in 
the SSFR-mass relation for low mass galaxies. Towards higher 
masses, the metallicity difference (between z~0.7 and z~0.1, eqns. 
7 and 8 of Savaglio et al. 2005) becomes smaller, in fact crossing 
at log(M*/Mo)Ri 10.3^, and so differential evolution in the M-Z 
relation cannot explain the upturn at log(A/,/A/0)^9 relative to 
log(Mt/M0)~lO galaxies. Furthermore, the evolution in average 
dust extinction between z~0.1 and z^l appears to be mild, even 
for the more massive galaxies considered here (e.g., Garn et al. 
2010). Thus, the expectation is that the empirical mass-dependent 
correction for [O II]-SFR should not change significantly by Z'^l. 
To properly test our empirical correction would require an analogue 
of our local test using B aimer decrement corrected Ha-SFRs at 
z~l, and such work is ongoing. 

5.2 Sample selection 

Another possible consideration is how star-forming galaxies are se- 
lected. It is now well-established that galaxies show a bimodal dis- 
tribution in colour out to at least z~l (Strateva et al. 2001; Baldry 
et al. 2004; Bell et al. 2004), comprising a 'blue cloud' dominated 
by star-forming galaxies, and a red-sequence dominated by passive 

^ This may not be physical and may be due to the use of a linear fit to 
the z~0.7 data, and a polynomial fit to the z~0.1 data. Nevertheless the 
metallicity difference at these masses is still small. 



galaxies with some minority of dust-reddened star-forming galax- 
ies. A similar bimodality is also seen locally in the distribution of 
SFR or SSFR (Brinchmann et al. 2004; Salim et al. 2007, McGee 
et al., MNRAS submitted), such that galaxies exhibit a relatively 
tight star-forming sequence (which may have scatter as low as 0.2 
dex (Salim et al. 2007) when plotted against stellar mass), and a 
broader population of galaxies with low (S)SFRs exhibiting little or 
no star-formation. However a tail exists down from the star-forming 
sequence such that it is necessary to go >1 dex lower in SFR to en- 
compass 95% of galaxies'^'' (fig. 17, Brinchmann et al. 2004). Thus, 
however galaxies have their SFRs measured, be it emission lines 
or SED-fitting (such as Salim et al. 2007, McGee et al., in prep), 
it is necessary to identify the star-forming sequence using either an 
overdensity selection in (S)SFR-mass space; or an approximate cut 
on (S)SFR; or, as we choose here, a colour cut to isolate the blue 
cloud. This rejects the tail of galaxies towards lower SSFRs and 
allows the bulk of the star-forming sequence to be isolated. This 
colour cut has the advantage when used with [O II]-SFRs that it re- 
jects objects for which the source of the emission is likely not star- 
formation, and it simplifies the comparison with other works where 
people have applied a colour selection. Throughout this work, we 
have applied a colour cut to exclude red galaxies (Fig. 1). This 
makes negligible difference to any of our results for dwarf galaxies 
since only one galaxy out of 199 in ROLES is red by this criterion, 
and the low mass end in the SDSS data is unaffected by the colour 
cut (c.f. left and right panels of Fig. 3). At the high mass end, the 
z~l data move slightly toward lower SSFRs (but still within the 
1(7 errors) when red galaxies are included in the computation of the 
mean (c.f. left and right panels of Fig. 3), suggesting that LINER 
activity and/or dust-reddened star-formation has little impact on the 
FORS2 sample. At z~0. 1, a significant lowering of the average 
SSFR in the SDSS Stripe82 data occurs when red galaxies are in- 
cluded in the mean. This is caused by the wide tail of low SSFR 
galaxies extending from the main star-forming sequence which can 
be seen as the difference between the two panels of Fig. 3. Exclud- 
ing red-sequence galaxies better traces the peak of the contours (i.e. 
the mode) of the star-forming sequence and is more directly compa- 
rable to the selection made at z~l, and so we favour this approach. 



6 CONCLUSIONS 

We have presented results from a spectroscopic sample of unprece- 
dented depth at z~l, probing star-formation via [O II] in galaxies 
down to log(M«/M0)~8.5. This is combined with an equally- 
deep comparison sample taken from Stripe 82 of the SDSS. We 
have, for the first time, measured the SSFR-mass relation and mass 
function of star-forming galaxies down to log(M,/M0)~8.5 at 
z~l using spectroscopy. 

The strength of the ROLES dataset is that it offers a highly- 
complete spectroscopic sample for such low mass galaxies. Mea- 
surements for these dwarf galaxies are particularly robust since they 
are numerous enough to provide useful statistics, and systematic 
uncertainties in SFRs due to extinction are negligible since they 
possess negligible dust, as confirmed by the lack of 24/im emission 
(Appendix B2.2). 

Dwarf galaxies are the building blocks of larger galaxies in the 
hierarchical formation scenario. The combination of ROLES and 



" where galaxies with spectroscopic indications of AGN activity are ex- 
cluded 
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Stripe82 data extends these observations at unprecedented depth 
to a wide redshift baseline, allowing evolution to be examined in 
detail for this important population of galaxies. 

The mass function for star-forming galaxies has remained re- 
markably constant over this interval, dropping by ~(6-35)% from 
z=1.0 to 0.1. This confirms that galaxies are leaving the star- 
forming sequence/blue cloud. 

The SSFR-mass relation decreases in normalisation towards 
the present day. Although at high masses this relationship evolves 
almost in parallel, the evolution is greater toward lower masses 
(log(M«/MQ);^9.5). The evolution of this relation can be ap- 
proximately modelled by the staged galaxy formation toy model 
of Noeske et al. (2007a) in which galaxies' SFRs decline expo- 
nentially with a timescale dependent on their baryonic mass (r oc 
M^^), and also the formation redshift depends on the galaxy bary- 
onic mass 0- + Zf oc M°'^). This low mass upturn can also be seen 
in SFRD-mass by looking at differences in the cumulative distribu- 
tions. Above log(M,/AfQ)>10, the z=0.1 and z=1.0 SFRDs are 
extremely similar, both containing ~70% of the cosmic SFRD at 
this epoch above this mass. This upturn is subtle, and at the limit 
of our data, but we have investigated and shown that it is unlikely 
to be due to incompleteness, or evolution in the mass-metallicity 
relation or dust content of galaxies. 

The Bower et al. (2006) version of the galform semi- 
analytic model of galaxy formation makes predictions for low mass 
galaxies which are in reasonable agreement with the measured 
SSFR-mass and SFRD-mass relations. At the lowest redshifts, the 
highest and lowest mass model galaxies exhibit too high a ratio 
of past-to-present SFR, possibly suggesting that star-formation has 
been quenched too early in the model. At z~l, the measurements 
are most robust at the low mass end, and this suggests that the GAL- 
FORM model does not reproduce the upturn seen in the observa- 
tional SSFR-mass relation at these low masses. One possible cause 
may be the too efficient termination of star formation in satellite 
galaxies ("strangulation"). 

Future work will extend these results to higher redshift, which 
is the interesting epoch for downsizing as hinted at by the prelimi- 
nary z=2 results, using the ROLES' technique with the new gener- 
ation of red-sensitive CCDs. Work is ongoing in parallel to inves- 
tigate the reUability of the z~l empirically-corrected [O II] -SFRs 
using NIR-MOS. 
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Table Al. Correction factors to convert stellar masses, M,, and SFRs, 
to the BG03 IMF. Numbers give the values to be added to quantities in 
the original IMF to obtain the corresponding value for a BG03 IMF. For 
example, logM*_BG03 = logM*_KroupaOi — 0.08. These conversions 
have some dependence on star-formation history and metaUicity. See text 
for details. 



IMF 


log(M*)corr 


log('I')corr 


Kroupa (2001) 


-0.08 


-0.18 


Salpeter(1955) 


+0.11 


-0.26 


Kennicutt (1983) 


+0.04 


-0.37 



APPENDIX A: IMF CONVERSIONS 

Different workers often assume different IMFs (each of which is in- 
voked to fix various problems with the adoption of an earlier IMF, 
such as Salpeter 1955) when computing stellar masses and SFRs. 
These derived quantities, in principle, may have different depen- 
dancies on the IMF since they are sensitive to stellar populations 
with different lifetimes. Although, in detail, assumptions other than 
the IMF may play a significant role in the quantities derived (partic- 
ularly for the stellar mass-fitting where the stellar population syn- 
thesis models used, the star-formation histories considered, and the 
choice of priors, etc. likely dominate the uncertainties, e.g., March- 
esini et al. 2009), one can attempt to correct for systematic offsets 
in a broad sense. Derivations from the literature for IMF transfor- 
mations (e.g.. Bell et al. 2003, SavagUo et al. 2005, Baldry et al. 
2008) are combined with conversions based on PEGASE.2 models 
as outlined in Gilbank et al. (2010a) to produce conversion factors 
from several of the most commonly-used IMFs (all the ones used 
in this paper: Salpeter 1955; Kroupa 2001 ; Kennicutt 1983) and the 
BG03 IMF. These are listed in Table Al. The SFR conversion fac- 
tors are based on the relative luminosity of the Ha line in scenarios 
with different IMFs, and are thus directly applicable to [O II] which 
is calibrated through the empirical relation with Ha. Another dis- 
tinction not always stated on the literature concerns which stellar 
objects are included in the estimate of mass. For example, ROLES 
and F0RS2 derived stellar masses at z~l (Paper 2) only include 
main sequence and giant stars and do not include the mass locked 
up in stellar remnants. Including or omitting remnants makes ~Q. 1 
dex difference (for the BG03 IMF), assuming solar metallciity and 
a population age of 10 Gyr, smaller than the random mass errors of 
0.2dex (Glazebrook et al. 2004). 

It is worth emphasising that the empirical mass-dependent cor- 
rection to [O II] derived by Gilbank et al. (2010a) is based on stel- 
lar masses and Snis assuming the Kroupa IMF. In order to apply 
eqn. 8 of Gilbank et al. (2010a), stellar masses and SFRs must first 
be converted to the Kroupa IMF and then back to the IMF of choice. 

APPENDIX B: COMPARISON WITH OTHER z~l 
[O II] SPECTROSCOPIC SURVEYS 

Bl Additional samples 

Additional z~l data are taken from DEEP2 spectroscopy of 
the AEGIS field (Noeske et al. 2007b, kindly provided by K. 
Noeske). These data are restricted to the region covered by MIPS 
24/Ltm imaging and where the MR photometry was deeper than 
K = 22. Stellar masses were fitted by Bundy et al. (2006). We 
have restricted the redshift range of this data to that of ROLES 
(0.88<z<1.15) and so the emission line estimates from DEEP2 all 
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come from [O II]. Unlike the above data, SFRs were not estimated 
purely from [O II] luminosity, but from the published data come 
from a combination of [O II]+24/im luminosity where 24/Ltm emis- 
sion is significantly detected, and from [O II] otherwise. So, some 
care is needed when comparing these different indicators. 24/im 
SFRs were computed by template-fitting (Le Floc'h et al. 2005; 
Chary & Elbaz 2001); and emission line SFRs were measured by 
Weiner et al. (2007). We have converted the assumptions used by 
the DEEP2 team in computing the purely-[0 II]-SFRs to be the 
same as ours (i.e., they assume [O II]/Ha=0.69 and an average 
extinction at Ha of 1.30 mag independent of mass/magnitude; we 
adopt a nominal [O II]/Ha=0.50 and 1 mag of extinction at Ha, 
to which the empirical mass-dependent correction of Gilbank et al. 
2010a can then be applied directly). However, we cannot separate 
out the contribution of the two components in the [O II]-|-24/im 
SFRs." Thus, the DEEP2 data consist of [O II]-SFRs using our 
mass-dependent correction (where no 24/xm detection exists) or an 
uncorrected [OII]+24/im SFR (where a significant 24/im detec- 
tion exists). The DEEP2 SFRs are only calculated for blue galaxies 
(Noeske et al. 2007b), since the emission line luminosity for red 
galaxies primarily comes from AGN/LINER activity rather than 
star-formation. All measurements of SFR and stellar mass are con- 
verted to our adopted BG03 IMF (see below). 

As a check of the other two z~l datasets, [O II] SFRs from 
the zCOSMOS survey (Maier et al. 2009, kindly provided by C. 
Maier) are used. These cover a slightly different redshift range 
(0.7 < z ^ 0.9) from that of the ROLES and DEEP2 data, and 
the effect of evolution over this range is examined in §B3. The 
Maier et al. (2009) SFRs were calculated using the empirical Mb- 
corrected formula of Moustakas et al. (2006). This is a better es- 
timate of the SFR than just using [O II] luminsoity, but still sig- 
nificantly underestimates the SFR at high stellar masses (Gilbank 
et al. 2010a). So, to correct for this and to measure SFR in a way 
consistent with the other samples, the Ms -corrected [O II] SFRs 
are multiplied by the ratio between the mass-dependent corrected 
[O II] SFR and the Ms-corrected SFRs at z~l from fig. 15 of Pa- 
per 2^^. These corrections are ~0.2-0.3 dex at log(Mt/M0)~ll. 

The above two datasets will next be compared with our data 
in the CDFS. To recap, all galaxies defined as 'star-forming' in 
these datasets, have a spectroscopic redshift, significant non-zero 
[O II] emission, and lie in the blue cloud. In the case of a subsample 
of the DEEP2 data, some also possess significant 24/im emission 
(but may or may not additionally exhibit [O II] emission). 

Bl.l 24 /dm data 

In order to fairly compare the combined [O II]-|-24/im SFRs from 
Noeske et al. (2007b) with SFR estimated solely from empirically- 
corrected [O II], we turn to 24/im observations of our CDFS sam- 
ple. The CDFS possesses deep 24/im MIPS observations from 
Spitzer (Dickinson et al., in prep.). MIPS 24/itm fluxes are taken 
from the FIREWORKS catalogue for objects in the spectroscopic 

Ideally, to correct [O 24/im SFRs, one would Uke to apply the 
empirical correction of Gilbank et al. (2010a) (which corrects for mass- 
dependent extinction and metallicity trends) to the [O II], and then divide 
this by the mass-dependent extinction (eqn. 9 of Gilbank et al. 2010a) to 
remove the dust correction, to which one may add the 24/im flux which 
accounts directly for the extinguished component of star-formation . 

This is actually the ratio of the SFBlDs but the data use exactly the same 
galaxies. In this way, the ratio gives a volume-weighted correction between 
the two SFR estimators. 



range of interest (0.88 < z ^ 1.15). One of the main disadvan- 
tages of MIPS photometry is that the PSF is relatively large (6") 
and so faint sources quickly become confused. To irunimise this 
effect, we begin by just considering high mass sources, which will 
generally have higher SFRs and hence should be the brighter 24/im 
sources. To achieve this, objects are selected from just the F0RS2 
sample. In order to attempt to quantify the impact of confusion 
(which may still be problematic for these brighter sources), ob- 
jects with significant IRAC band 1 flux (3.6/im, hereafter, [3.6]) 
are selected froin the FIREWORKS catalogue. The [3.6] flux is of- 
ten used as a prior when attempting to deblend the 24/xm images, 
since tj^jically every 24/im source is detected in [3.6] and the PSF 
of the latter is much better than the former (1.6"). The [3.6] source 
list is cross-correlated with the 24/tm list, and objects with more 
than one match with 6" are flagged as potentially blended. The 
6" radius means that two lAjim sources separated by this amount 
will touch at their half-maxima. This should only be regarded as 
a crude estimate of confusion since the FIREWORKS catalogue 
is if -band selected, and so a number of fainter [3.6] sources do 
not appear even in the deep if-band imaging, and hence some po- 
tential blends may not be identified. In addition, the presence of a 
second [3.6] source within the matching radius does not guarantee 
it will be a significant source of 24/im flux, and so cases flagged 
as 'blended' may in fact suffer no contamination from a neigh- 
bour. Nevertheless, this method provides an approximate division 
between blended/unblended sources. 



B2 Comparison of [O II] with 24/im SFRs 

B2. 1 High mass galaxies 

To understand any differences on our results from comparing 
our empirically-corrected [O II]-SFRs with [O II]+24/im SFRs 
from Noeske et al. (2007b), the 24/im data are considered for 
the bright/high mass (and likely least-confused) FORS2 galaxies 
in CDFS, as described in §B1.1. Keeping the previously men- 
tioned caveats regarding source confusion in mind. Fig. B 1 shows 
the comparison between [O II], emp.corr -SFR and [O II]-t-24/im 
SFRs, where the [O II] has not been empirically corrected, nor cor- 
rected for any extinction, before adding to the IR luminosity. For 
consistency with the DEEP2 method, the assumed [O II]/Hq: ratio 
used here is 0.69 (rather than 0.5) and no correction for extinc- 
tion is made (since the light from star-formation being reprocessed 
by dust should now be measured directly by the 24/im flux. The 
[O II] SFR using the DEEP2 assumptions is denoted SFR([0 II] d) 
to be explicit. To convert the observed 24/im luminosity to total 
IR luminosity, the templates of Chary & Elbaz (2001)^^ are used 
(as used by Noeske et al. 2007b). The total IR luminosity is then 
converted to SFR using the prescription of Bell et al. (2005) con- 
verted to our IMF. Filled circles show galaxies which are not con- 
sidered confused, using the above (strict) criteria, and open circles 
show those which would be considered confused. Crosses indicate 
galaxies which do not have significant 24/im detections, and so the 
[O II]+24/tm SFR comes solely from [O II]. As can be seen, the 
majority of these high mass objects show significant 24/im emis- 
sion. The green square indicates a source which may be an AGN 
based on its MIR colours (see Appendix C), so the 24/im flux may 
be contaminated by non-thermal emission. The blue diamonds de- 
note objects which may contain an AGN based on X-ray emission 



using routines from: http://david.elbaz3.free.fr/astro.codes/chary_elbaz.html 
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Figure Bl. Comparison between empirically-conected [O II]-SFRs and 
[O II]+24/^m SFRs, as used by DEEP2/ AEGIS. Circles show all measure- 
ments with significant 24/^m detections; crosses indicate S/N(24/im)<3. 
Filled circles show which galaxies are isolated by our automated criterion 
based on checking for [3.6]-detected neighbours in the _ff -selected FIRE- 
WORKS catalogue. Blue diamonds denote X-ray detections, which sug- 
gests that the [O II] and lA/im fluxes are contaminated by an AGN, and the 
green diamond indicates the possible presence of an AGN based on MIR 
colours (see Appendix C). Note that the isolation criterion is likely too strict 
(as discussed in the text), and the majority of open circles are likely usable 
measurements, as supported by the fact that most open and filled circles 
follow the same average trend. 

(Appendix C). The dashed line shows the best-fit relation, which 
is offset such that the SFR([OII]d + 24/im) estimate is a factor 
of 2.2 higher than SFR([OII], cmp.corr.). Using instead the local 
conversion of Rieke et al. (2009) results in an even greater offset 
of the SFR including 24/im data. Without another independent es- 
timate of the SFR, we cannot assess whether this offset is due to 
the empirical [O II] estimate underestimating the total SFR, or the 
[O II]+24/^m measurement overestimating it. To first order, we are 
primarily interested in the comparison between [O II] -based SFRs 
at z=0.1 and z=1.0. The relative difference between the purely 
[O II] -based measurements and those including 24/im data at z=l 
will be considered below. 

B2.2 Stacking analysis for low mass galaxies 

Following the same method as for the higher mass FORS2 sub- 
sample above, isolated, lower-mass galaxies from ROLES are se- 
lected. Seven ROLES galaxies classed as isolated are individu- 
ally detected in 24/im emission, and they follow the same rela- 
tion, within the broad scatter, as the higher mass galaxies above. 
Considering only isolated galaxies without individual 24/im detec- 
tions results in a sample of 32 objects. The 3a limit of the 24/im 
photometry in CDFS is 11/iJy (Wuyts et al. 2008), which corre- 
sponds to a SFR of O.TMq yr~^ at z=l. The fact that the major- 
ity of isolated ROLES galaxies are undetected down to this limit 
is reassuring that these low mass objects possess little dust, and 
hence the [O II] line should be tracing reliably their SFRs. In or- 
der to push to deeper limits, cut-outs are extracted from the 24/im 
image (obtained from the GOODS-S v0.30 release, Dickinson et 



al., in prep.) around the K-bznd position of the undetected objects 
and the mean of the stack is found. The total flux is measured in 
a series of apertures ranging from 3" to 6". No significant detec- 
tion was made down to the limit of the image, which should be 
2 (=Il/\/32) fJ-Jy- Errors were calculated from 100 bootstrap re- 
samplings of the 32 images used in the stack. The formal limit 
comes out to be —0.17 ± 0.24/iJy. This is many times lower than 
our [O II]cmp.corr limit (SFR>O.3M0 yr~^), implying that the vast 
majority of star formation is unextinguished in these mass systems. 

B3 The SSFR-mass relation 

As in Fig. 3, SFRs and stellar masses from the z~O.I and z~l sam- 
ples are combined to study the evolution of SSFR-mass in Fig. B2. 
The two panels show the different z~l data. In both panels; the 
SDSS (z~0. 1) galaxies are shown as contours for clarity. Filled red 
circles show the mean in bins of stellar mass. For the left panel, 
the effect of including all galaxies rather than just those which be- 
long to the blue cloud is shown by the dotted contours with the 
mean relation as open red circles. This would give a lower aver- 
age SSFR at higher masses, since higher mass galaxies possess a 
long tail toward lower SSFRs. The z~l data are shown as smaller 
filled circles, colour coded by dataset: black points are ROLES' 
galaxies, blue points are blue, emission line only galaxies from 
DEEP2 and orange circles are DEEP2 [OII]+24/im SFR mea- 
surements. Larger black circles show the median in bins of stel- 
lar mass for the z~l data. Dashed lines indicate SFRs of 0.1 Mq 
yr^^, 1 M0 yr~^, and 10 Mq yr^^. The black dotted line indicates 
a power law fit to the z=l data which is given by \og{SSFR) = 
—0.138 log(A/*) — 7.95. As can be seen, this fit is consistent with 
the bulk of the points but does not capture the turnover at the high- 
mass end, log(Af,/A/0)j>lO.5. [The higher mass points are better 
fit with log(5'5'Fii) = ^.5771og(M*) - 3.31.] 

The green curves show the best fit tau model (as described 
in §4.1) for these data, with parameters as annotated on the plot. 
The best fit using the DEEP2+R0LES data shows a much flatter 
dependence, almost zero, of z/ on Alt (P = 0.06) than from the 
F0RS2+R0LES data (§4. 1, and repeated in right panel of Fig. B2). 
The data points between lO<log(M,/M0)<l 1 were not consid- 
ered in Noeske et al. (2007a) due to their <95% completeness but 
including or removing these makes little difference to the fit. 

The right hand panel of Fig. B2 shows the same information, 
but the DEEP2 data have been replaced by the data from the public 
F0RS2 spectroscopy in CDFS (smaller red filled circles). Recall 
that, for all these data, only blue cloud galaxies are used, and al- 
though individual data points are shown for each galaxy, the mean 
is calculated weighting each galaxy by a completeness and Vmax 
weight. For the F0RS2 data, the mean relation is not significantly 
lowered by including all galaxies with [O II] emission, rather than 
just blue cloud galaxies with [O II] emission. 

Clearly there is significant disagreement at the high mass end 
between the two different surveys used. To check that the differ- 
ence between our R0LES+F0RS2 dataset and the DEEP2/AEGIS 
data is not due to cosmic variance from our relatively small vol- 
ume, we compare our measurements with the larger zCOSMOS 
survey. zCOSMOS SFRs from [Oil] (Maier et al. 2009), cor- 
rected to use our empirical mass-dependent correction, are taken. 
These data (plotted as thick open blue diamonds with error bars in 
the right panel of Fig. 3) cover a slightly different redshift range 
(0.7 < z < 0.9) from that of the ROLES and DEEP2 data. Esti- 
mating the size of this evolution by measuring the evolution from 
the DEEP2 data cut to the same redshift bins gives an increase of 
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Figure B2. log(SSFR) versus log(stellar mass) for galaxies in our samples. Small, filled points show data at 0.88 < z ^ 1.15: black points are from ROLES, 
blue points are emission-line only SSFRs from DEEP2 and orange points are [O II]+24/im SSFRs from DEEP2. Contours denote SDSS data at z~0.1. Larger 
filled circles show mean SSFRs in bins of stellar mass for the data in the two redshift ranges: black points at z~ 1 and red points at z~0. 1 . The dotted contours 
and open red circles show how the SDSS data would change if all galaxies with significant H« flux were included, rather than considering only blue cloud 
members, due to the long tail toward lower SSFRs at higher masses. Dashed lines denote the three different SFRs, as annotated. In the right panel, the higher 
mass data (smaller red filled points) are taken from ESO public spectroscopy rather than DEEP2. Thick blue diamonds with error bars show results taken from 
zCOSMOS (0.7 < z 0.9); thinner blue diamonds show these results evolutionary-corrected to 0.88 < z 1.15 as described in the text. Dotted black 
lines show best fit power laws to the z=l data, as described in the text. Green lines show best-fit tau models, discussed in §4.1, with parameters as indicated 
on plot. See text for details. 



0.15-0.25 dex between 0.7 < z ^ 0.9 and 0.88 < z < 1.15 in the 
two mass bins. This moves the open diamonds in Fig. 3 into closer 
agreement with the CDFS data (thinner blue diamonds). Thus, we 
are reassured that our high mass points in CDFS are not severely 
affected by cosmic variance. 

At the high mass end, say log(Mt/M0)~lO.6, the aver- 
age SSFR measured in DEEP2 is 0.5 dex higher than that mea- 
sured in CDFS or zCOSMOS. Part of this discrepancy may be 
due to the use of different SFR indicators. Fig. B 1 shows that us- 
ing SFR([0 Il]D+24^m) instead of SFR([0 lI]),cmp.corr. leads to 
a factor of 2.2 (0.35dex) higher estimate of the SFR. Kelson & 
Holden (2010) have suggested that using observed 24fim at z~l 
to estimate SFRs may overestimate the ongoing SFR by factors of 
~1.5-6 since dust-encircled thermally pulsating-asymptotic giant 
branch (TP-AGB) stars with ages of 0.2-1. 5Gyr may constitute a 
significant source of contamination to flux interpreted as coming 
from star-formation. Early results from the Herschel Space Obser- 
vatory show that there is considerable scatter between TIR lumi- 
nosities measured directly from near the IR peak of the SED (at 
lOO^tm and 160/im) and those extrapolated from observed 24/im 
luminosities, with a systematic bias which is likely a function 
of redshift (Rodighiero et al. 2010). Nordon et al. (2010) find 
that Spitzer 24/^m observed SFRs overestimate /ferac/jeZ-measured 
SFRs by factors of --^4-7.5 at slightly higher redshift (1.5<z<2.5). 
The details of these likely depend on the SED-fitting method to 
both datasets, and it is not trivial to see how easily these re- 
sults transfer to the data considered here. However, even consid- 
ering just the DEEP2 points using only [O II] SFRs (generally 
lower mass, log(Af»/Af0)~lO) shows a significantly higher SFR 
in the[0 11] data versus the CDFS data, and so the difference cannot 
be attributed entirely to the inclusion of the 24^m data. 

One main difference between the DEEP2 [O II] data and the 
other two datasets is that DEEP2 is not flux calibrated (as the 
other samples are), but line strengths are calculated using equiv- 
alent widths and broad band photometry to calibrate the contin- 



uum. Another difference is that DEEP2 is very biased towards blue- 
selection relative to the other surveys (CDFS is z-band selected, 
zCOSMOS is J-band selected). Although the A'-band follow-up 
photometry helps with the mass selection of the AEG1S-DEEP2 
data, the initial DEEP2 spectroscopy is still i?-band selected (which 
is close to rest-frame B at z~l, and thus more closely resembles a 
SFR-seleted sample rather than mass-selected). Thus, galaxies with 
low SFRs might possibly be missed, leading to an absence of low 
SSFR objects. In practice, it may be that a subtle interplay of fac- 
tors are responsible for the disagreement of the DEEP2 data with 
the other results. Indeed, Cowie & Barger (2008) make a similar 
suggestion when arguing why they find a much longer tail toward 
lower SSFRs compared with the DEEP2 results. 



APPENDIX C: AGN 

In order to test for contamination of the [O II] flux or 24/im flux 
by AGN in our sample, the technique of Stem et al. (2005) is used 
to examine the MIR colours for the signature of an AGN. [5.8]- 
[8.0] vs [3.6]-[4.5] colours for all galaxies in both the ROLES and 
F0RS2 subsamples of CDFS are shown in Fig. CI. Objects with 
greater than 5cr detections in each of the four bands are indicated 
as filled circles with error bars. Only these significant detections 
are considered when testing for the presence of an AGN, as reli- 
able colours cannot be measured for the other objects. The dotted 
line indicates the region within which AGN should lie according to 
the criteria of Stern et al. (2005) for galaxies over the redshift range 
Q < z < 2. We refine these criteria slightly considering the bound- 
ary around AGN-like objects in the redshift range 0.8 <; 2 1.2 
by reading values off figs. 2 & 3 of Stem et al. (2005) which bracket 
the majority of the AGN population. This results in the green lines 
plotted in Fig. CI. As can be seen, only a small minority of galaxies 
(with small photometric errors) lie within this region . AGN follow- 
ing these criteria are indicated by green squares. Only 15 galaxies 
out of our total sample of 3 1 1 meet this definition, and only two of 
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Figure CI. MIR colour-colour diagram to examine for the presence of 
AGN. Points are taken from both ROLES and F0RS2 spectroscopy in 
CDFS. Objects with > 5cr detections in all four filters are shown as filled 
circles with error bars. Red open circles indicate > 3(t 24 /im detections. 
The dotted line indicates the region populated by AGN according to the gen- 
eral < 2; < 2 criterion of Stem et al. (2005). The solid green box shows 
a refinement of this designed to target our redshift range. Blue diamonds 
indicate X-ray detected galaxies. 

these are drawn from the FORS2 sample, plotted in Fig. Bl. It is 
also instructive to examine which of these objects are X-ray lumi- 
nous and thus classified as AGN this way. Four such X-ray detec- 
tions are present in the combined R0LES+F0RS2 sample (three 
from F0RS2; one from ROLES). Of these four, one is undetected 
in all four IRAC bands (indicated as the point at 0,0 in the plot), 
and one is only significantly detected at [5.8]. The remaining two 
are shown in the plot and lie just outside our z~l AGN criteria: one 
is within la of the edge, and the other just over 2a away. This sug- 
gests that the MIR colour is not an unreasonable method of select- 
ing AGN. More importantly, this test suggests that contamination 
of our SFR estimates by AGN is likely to be minimal, and we are 
justified in neglecting it in our analysis. 
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